Trapped 6 Li : A high T c superfiuid ? 
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We consider the effect of the indirect interaction due to the exchange of density fluctuations 
on the critical temperature of superfiuid 6 Li . We obtain the strong coupling equation giving this 
critical temperature. This equation is solved approximately by retaining the same set of diagrams as 
in the paramagnon model. We show that, near the instability threshold, the attractive interaction 
due to density fluctuations gives rise to a strong increase in the critical temperature, providing a 
clear signature of the existence of fluctuation induced interactions. 



PACS numbers : 32.80.Pj, 67.90.+Z, 74.20.Fg 

The recent quite impressive progress in obtaining ul- 
tracold atomic gases has open the way to the discov- 
ery of a large number of new superfluids. Already most 
alkali Bose gases have been shown to undergo Bose Ein- 
stein condensation. Even if superfluidity has not been yet 
demonstrated explicitely in experiments, they are firmly 
believed to be superfluids because they have already been 
seen to display phase coherence. Moreover the search for 
the transition of trapped Fermi gases toward a BCS type 
superfiuid is now actively considered 0] , as the possibility 
of experimental observation is quite realistic. 

A high critical temperature will be obtained for a 
strong attractive effective interaction between fermions. 
Since at low temperature scattering is essentially s-wave, 
a large negative scattering length is most favorable. 
Therefore spin polarized 6 Li appears as a strong can- 
didate since its scattering length is found experimentally 
to be very large a = - 1140 A . In this case pair- 
ing would occur between two different hyperfine states 
0, which play the same role as spin states in standard 
BCS theory. This is the situation we will consider and we 
will assume the most favorable case where the number of 
atoms is the same in the two hyperfine states. It is natu- 
rally important to assess as well as possible the value of 
the critical temperature. Experiments are presently per- 
formed in magnetic atomic traps leading to a harmonic 
potential and calculations should be done for this geome- 
try |H . Or one can also consider making use of an optical 
trap. Here we will restrict ourselves to the simpler sit- 
uation of fermions confined in a box with total density 
n = k 3 F /3n 2 . 

The high density regime is of particular interest since 
it corresponds clearly to higher critical temperature and 
the superfiuid will be more accessible experimentally in 
this regime. However there is also a deep theoretical 
interest to investigate this domain. Indeed the high den- 
sity regime is bounded by an instability which occurs 
for a coupling constant A = 2kp\a\/TT > 1 . Beyond 
this limit the compressibility becomes negative because 
of the strong effective attraction between atoms in dif- 
ferent hyperfine states. In the vicinity of this instability 



A < 1 , the compressibility becomes high and density fluc- 
tuations occur easily. This leads to the possibility of an 
attractive interaction between fermions through the ex- 
change of density fluctuations, in a way completely anal- 
ogous to the phonon exchange mechanism of standard 
superconductivity. Qualitatively we expect this mecha- 
nism to add up to the direct attractive interaction and 
to lead to an increase in the critical temperature. 

Independently of this increase, this situation is ex- 
tremely interesting because it is quite analogous to what 
is believed to happen in many other condensed matter 
systems. High T c superconductors and the newly dis- 
covered low temperature superconductor Sr2RuOi are 
well-known examples. However the best case is prob- 
ably liquid 3 He Q which is not so far from having a 
ferromagnetic instability. It has been proposed that 
strong spin fluctuations, the so-called paramagnons, ex- 
ist in this liquid and play an important role in the 
physics. In particular the attractive interaction, leading 
to the Cooper pairs formed in the superfiuid, has been 
attributed to paramagnon exchange. A strong qualita- 
tive support for this picture is the existence of the A 
phase of superfiuid 3 He at high pressure, which has been 
explained by feedback effects of its specific structure 
on the paramagnon propagator. However it is not clear 
at all that the paramagnon model can provide a quan- 
titative description of the properties of the liquid, and 
in particular can account for the observed values of the 
critical temperature || . 

In this respect the situation in 3 He is difficult, since one 
does not have a precise knowledge of the instantaneous 
part of the pairing interaction. Moreover the parameter / 
involved in the paramagnon description varies in a quite 
restricted range in the vicinity of the instability limit I 
= 1, when the pressure of the liquid is varied over the 
full range available in the phase diagram. By contrast 
the situation potentially offered by superfiuid 6 Li gas is 
much more agreable : the instantaneous interaction can 
be directly linked to the diffusion length, which is fairly 
precisely known. Moreover the possibility of varying to 
a large extent the density allows to change the coupling 
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constant at will. This offers a stringent coherence check 
of the theory since experiment will be able to verify the 
general qualitative behaviour of the model. In this way 
we can hope to have a definite answer to the question : 
is the paramagnon model a proper description or not ? 

The problem of the critical temperature for a BCS su- 
perfluid in a dilute Fermi gas has been investigated by 
Gorkov and Melik-Barkhudarov ||, following the work 
jlO| of Galitskii. It has been recently considered in the 
more general context of dilute atomic gases Q by Stoof 
et al. who found a typical value of 40 nK for a density 
n = 10 12 cm~ 3 . This temperature is quite within reach 
experimentally. Our purpose is to extend these works in 
the high density regime. Naturally our result will reduce 
to the proper one in the dilute limit. 

The critical temperature is obtained ^ by writing that 
the vertex part in the normal state diverges, which is ex- 
pected to occur first for zero total momentum and en- 
ergy of the pair. The corresponding vertex part T p y 
is related to the irreducible vertex r p p / in the particle- 
particle channel by : 

r P:P ' = r PiP ' — f Ptk D k T kiP ' (1) 
fe 

where J2 k 1S f° r J" c ^ c ' We have set D k = 

G(k)G(—k) where G(k) is the full Green's function and 
k = (k, u> n ) is a four-momentum. The summation runs 
over the wavector k and the Matsubara frequency uj n — 
(2n+l)nT. 

We split the irreducible vertex into the bare interac- 
tion £7 PjP ' and all the contributions T* p , which are higher 
order in the interaction : f p y = C/ p , P ' + T* , . Then, 
following Galitskii (To), we eliminate the interaction U 
in favor of the vertex in the dilute limit corresponding 
physically to two atoms scattering in vacuum. In this 
limit, f reduces to U, and G(k) becomes the free partic- 
ule Green's function (ito n — t k )~ l with e k = k 2 /2m. In 
contrast to Ref. Q , we have taken for convenience the 
chemical potential equal to zero in this limit. Let us call 
Tp , the vertex in this limit. From Eq.(l) it satisfies : 

where D k = (uj 2 + e 2 .)^ 1 ■ Since we will deal with 
fairly small temperature, we can take the T — * limit 
where r p , reduces to the vertex T p , for two scatter- 
ing atoms evaluated at zero energy. This vertex can be 
explicitely expressed pJj ] in terms of the scattering ampli- 
tude /(p, p') corresponding to the scattering potential U. 
Since the atomic potential has a very short range com- 
pared to all other lengths involved in the problem, the 
typical wavevector for a change in £/(p,p') is very large 
compared to the wavevectors we have to deal with and 
similarly for T p , . Hence we can take in our problem 



T p , equal to its p = p' = limit, which is given |Tl| by 
r° = Ana/m in terms of the scattering length. This can 
be checked explicitely in the case of a separable potential 
or for a pseudopotential JL2| . 

Since T p y diverges at T c , the first term in the r.h.s. 
of Eq.(l) is negligible. In the resulting integral equation, 
p' appears as a free parameter which can be omitted. 
Writing T p instead of T p/p i , we obtain that T — T c when: 

r P = - r EP P ,k+r; fc pA (3) 

is satisfied. Note that the effective interaction T* k will 
be frequency dependent. Therefore we have also to retain 
the frequency dependence of Tk ■ In other words the order 
parameter in the superfluid phase will have a frequency 
dependence, which corresponds to a strong coupling sit- 
uation. We eliminate U from Eq.(3) by premultiplying it 
by <5p', P — TT p , p D p and summing over p. Making use of 
Eq.(2) this leads to: 

r p = -r°Tj2[Dk-D° k }r k 

k 

-TY^K,k-r Tj2 D k> r t<,k}D k r k (4) 

fc k' 

Since we expect the dependence T p to be fairly slow, 
over a typical scale fejr, with respect to the wavevector p, 
we will neglect it in the following. For coherence T* k in 
the r.h.s. of Eq.(4) will be evaluated for |p| = kp. Sim- 
ilarly we can evaluate the self-energy S(fc) at the Fermi 
surface. In terms of the renormalization function Z n de- 
fined by Yi{kp,iui n ) — E(fejr,0) = iuj n {\ — Z n ) this leads 
to D k = (w*Z* + with Cfc = e fe - E F . This pro- 

cedure is the standard one for classical superconductors, 
fully justified by the existence of a small energy scale 
which is the phonon energy. Here to the contrary we 
have the single energy scale Ep , and the quality of this 
set of approximations is less obvious. Nevertheless, since 
one does not expect any drastic effect to arise from this 
k dependence and that this procedure is consistent with 
the paramagnons approach, we will work with this ap- 
proximate scheme. 

Once r p has no p dependence, integrations over mo- 
mentum can be performed in Eq.(4). In the first term, 
we write D k - D° k = [D k - D\] + [D\ - D° k ] where 
D l = ( w « + ^ly 1 - Tnc contribution from D\ - D° k can 
be integrated exactly. For D k — D\ we take into account, 
consistently with our above approximation, that uj n is 
small. In this way we obtain : 

/ dk [D k - D*] = [Cf - C7f ] (5) 

where Nj — mkp/2ir 2 is the density of states at the 
Fermi surface. We have set C* f = 1 - Z" 1 and = 

[ V /f + (l+ W 2) 1 /2_ y^li/yS With W n = UJn/Ep. OllC 
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can check that the term gives in Eq.(4) a con- 
tribution r°iV/ln(8e c - 2 .EF/7rT) (C is the Euler con- 
stant). If only this term is retained one obtains T c / Ep = 
8e c - 2 /7rexp(-l/A) used in 0. The term will give 
a decrease of the critical temperature due to mass renor- 
malization and lifetime effects. 

We turn now to the irreducible vertex T k k , . As we 
have indicated we are mostly interested in the contribu- 
tion of density fluctuations to this vertex, and here we 
will handle it by retaining the same set of diagrams as 
in paramagnon theory. Actually since the attractive in- 
teraction acts only between atoms with different hyper- 
fine states, we are exactly in the same situation as for 
paramagnons where interaction takes place only between 
different spins. The only qualitative difference is that 
the interaction is repulsive in paramagnon theory, leading 
to a positive dimensionless coupling constant I = Nfl, 
while the attractive interaction between 6 Li correspond 
to a negative I. Another important difference is that it 
would be inaccurate to retain only the bare interaction for 
all the elementary vertices in the paramagnon diagrams. 
Indeed we know that we have a large scattering length. 
This is obtained quantitatively by summing up the lad- 
der diagrams for two scattering atoms, as it is clear from 
Eq.(2). Obviously we have to do a similar summation 
in the paramagnon diagrams, otherwise we will miss the 
dominant contribution. More precisely we would need to 
know the irreducible vertex in the particle-hole channel. 
We will assume that the dominant contribution to this 
vertex is given by the sum of the ladder diagrams. Also 
we will not attempt to take into account its energy de- 
pendence and consider that this vertex is the same as for 
two atoms in vacuum. With these hypotheses we are led 
to take the interaction / of paramagnon theory equal to 
r°. This gives us I = N f T° = - A. 

Now the sum of the paramagnons diagrams ( includ- 
ing ladder and bubble diagrams ) gives j|Jl^l T k k , = 
N f V eff (k — k') with N f V cff {k) = J 2 xo/(l - Ixo) + 
7 3 Xo/(l — I 2 Xo) an d Xo(k) is the dimensionless elemen- 
tary bubble H . When the self-energy is evaluated with 
the corresponding set of diagrams, one finds M that it 
is given by the same expression as for an effective in- 
teraction N f V z (k) = Pxl/0- - Ixo) + / 2 Xo/(l - I 2 xl) 
In agreement with our above approximation, we evaluate 
the self- energy for wavevector kp and for small energies. 
This leads to: 

n 

(27i + l)(Z n - 1) = Vz(0) +2j2Vz{2npT) (6) 

P =i 

with V z (u p ) = (1/24) loQdqNfVziq^p). This V z (u p ) 
corresponds to the average of the interaction for scatter- 
ing of an atom at the Fermi surface from k to k' with 
wavevector transfer q = 2kp sin(0/2) where is the angle 
between k and k'. This angular average has to be per- 
formed numerically, in order to obtain Vz ■ As expected 



Z n > 1. The maximum is at zero energy, with a fairly 
long tail at high energy. Naturally Z n increases when 
|7| — > 1 and it diverges at zero energy in this limit. 

In the term F* k DkTk in Eq.(4), the only strong de- 
pendence on k, for low energies, comes from which 
forces k w hp- So we integrate over as above. Setting 
r„ = r(w n ), this leads us to: 

where as above V c &{uj p ) = — (l/2fc|,) j Q qdqNfV e g(q,w p ). 
Finally the contribution from the last term in Eq. (4) with 
the double summation is somewhat more complicated to 
evaluate. Nevertheless after a double integration, per- 
formed numerically without problems, which provides a 
function V c (u n ), it can be written as: 

t 2 £ D° k >n>, k D k r k = N f J2 f^^KO (8) 

Finally Eq.(4) can be cast as: 
r n = 2. 1 \^C m Y m + 7-=—V e g(uj n — uj m )T m (9) 

— Wm Y^raY^ra 
m 1 1 11 

where we have set C n = C r f - - C% with C% = 
V c (u> n )/Z n . Hence T c is the highest temperature for 
which the matrix corresponding to the r.h.s. of Eq.(9) 
has an eigenvalue equal to 1. Eq.(9) is very similar to 
the one obtained from Eliashberg equations for a strongly 
coupled superconductor. The most noticeable difference 
is the term produced by the interaction through fluc- 
tuation exchange when one replaces the bare potential in 
terms of the scattering length. We note also that, in the 
dilute limit, we recover the result of Gorkov and Melik- 
Barkhudarov M. 

Let us now consider the physical effect of V e g. It 
is more convenient to decompose it {j],|l3| as NfV e s — 
(3/2)J 2 xo/(l-/Xo)-(l/2)/ 2 Xo/(l+/Xo). The first term 
is due to spin fluctuations and the second one to density 
fluctuations. As pointed out by Berk and Schrieffer |J 
the first term gives always a repulsion which increases 
strongly in the vicinity of the ferromagnetic instability 
I — > 1. However in the range I < which is of interest 
to us, the density fluctuations attraction can take over 
the spin fluctuations repulsion. This is what happens in 
the limit I — > — 1 where the gas becomes unstable. Then 
in complete analogy with the paramagnon case, the ef- 
fective interaction at the Fermi surface V e f{(uj n ) diverges 
for zero frequency. For nonzero frequency it remains fi- 
nite. In particular at high energy where xo S e ts small, 
spin fluctuations dominate, NtV e g « -f 2 Xo an d the over- 
all interaction is repulsive. So for I near -1, the effective 
interaction is attractive at low frequency and has a re- 
pulsive high energy tail. But the attractive part exists 
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only for I < —0.6. For lower I fluctuations are repulsive 
for all energies. 

The divergence of V e R for J — ► — 1 may raise hope that 
in this limit T c is going to be very large. However it is 
known, in the case of strongly coupled superconductor 
as well as for [|| 3 He in paramagnon theory, that this in- 
crease is countered by the corresponding increase of Z n , 
physically due to mass renormalization and lifetime ef- 
fects. The net result of these opposite effects has to be 
found numerically, so we turn to our numerical results 
from Eq.(9) for T c . They are given in Fig.l. 

Surprisingly, instead of a regular rise of T c /Ep as a 
function of A = —I, we find two regimes. Up to A « 0.4 
we have a regular and strong increase (note that for 
A = 0.4 our result is markedly below the extrapolation of 
the result of Ref. §, which would give T C /E F = 0.023). 
We have then a saturation and even a slight decrease up 
to A w 0.6. We attribute this effect to the increase of Z n 
which overcompensates the increase with A of the direct 
attractive interaction. Then, starting for A rs 0.6 we ob- 
tain another strong rise of T c /Ep. This second regime is 
clearly due to the increasing contribution of the attrac- 
tive interaction due density fluctuations. T c /Ep grows 
up to a maximum 0.025 found for A ~ 0.98 (correspond- 
ing to T c — 190 nK). We stress that, if observed, this crit- 
ical temperature would be the highest (relative to Ep) 
among BCS superfluids, since for standard superconduc- 
tors as well as superfluid 3 He this ratio is of order 10~ 3 
whereas for high T c superconductors (if they are BCS) it 
reaches at best 10~ 2 . Then, when A increases further, 
T c decreases gently in the vicinity of A = 1 . This effect is 
clearly due to the increase of Z n . However the maximum 
of T c is obtained for A so close to 1 that this decrease is 
probably unobservable experimentally. 

The most interesting feature of these results is that the 
existence of an indirect attractive interaction due to den- 
sity fluctuations exchange gives a qualitative signature in 
T C (A). Although we may wonder about the quantitative 
validity of our results, due to the approximations we have 
made, it is reasonable to believe that the qualitative rise 
in T C (A) will survive. Its observation would be a strong 
indication of the importance of fluctuation exchange in 
the effective interaction in 6 Li. And indirectly it would 
also bring some support to the existence of similar mech- 
anisms in other BCS superfluids. Finally in this regime 
we expect, in the superfluid phase, deviations from stan- 
dard weak coupling BCS theory and feedback effects. In 



conclusion we have shown that, for 6 Li near the instabil- 
ity threshold, indirect attractive interaction due to den- 
sity fluctuations exchange can lead to rather high critical 
temperature, with a clear signature for its dependence as 
a function of the coupling constant. The observation of 
this effect would be very interesting, as a clear example 
of collective mode induced superfluidity. 
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FIG. 1. T c /Ef as a function of the coupling constant A. 
The dots correspond to actual calculations, the line is just a 
guide for the eye. 
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